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ABSTRACT 


High  pressure  shock  waves  were  studied  in  porous  copper  samples. 
These  waves  were  induced  by  high  velocity  impact  of  flat  plates. 
Flat  plates  of  copper  and  a  tungsten  alloy  were  projected  at 
velocities  up  to  8  mra/usec  using  a  two-stage  light  gas  gun. 
Samples  of  porous  copper  with  distention  ratios  of  3.08, 

4.46  and  8.74  (32.4%,  22.4%  and  11.4%  of  solid  density)  were 
studied  through  shock  transit  measurements.  Analysis  of  the 
experiments  reveaxs  that  a  complex  (two-wave)  shock  structure 
exists  in  the  open  celled  materials  impacted  at  conditions 
that  would  result  in  shock  stresses  up  to  4  Mb  in  solid  tar¬ 
gets.  The  nature  of  the  double  wave  precludes  the  use  of  the 
Rankine-Hugoniot  equations  which  assume  that  thermodynamic 
equilibrium  is  attained  behind  a  discontinuous  shock. 
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SECTION  I 

INTRODUCTION 


This  report  describes  work  performed  by  the  Materiel s  and 
Structures  Laboratory  of  Manufacturing  Development, .General 
Motors  Corporation,  under  Contract  No.  DAAD05-71-C-0403 
monitored  by  the  U.S.  Array ,  Aberdeen  Research  and  Development 
Center,  Aberdeen  Proving  Grounds,  Maryland.  The  purpose  of 
this  contract  was  to  determine  experimentally  the  IJugoniot 
equation  of  state  of  highly  porous  copper. 

There  are  two  areas  of  interest  surrounding  the  shock  compres¬ 
sion  of  porous  materials.  As  absorbers  of  mechanical  energy 
porous  materials  combine  the  desirable  qualities  of  good 
mechanical  dissipation  through  large  strains  at  low  stresses 
and  light  weight.  Secondly,  insights  into  the  physics  of 
materials  at  high  pressures  and  temperatures  may  be  obtained 
through  the  study  of  shock  compression  of  porous  solids . 

Many  studies  of  the  dynamic  response  of  porous  metals  to  sud¬ 
denly  applied  loads  have  been  performed  in  the  low  stress 
regime  up  to  about  50  kb.  Additionally,  a  good  deal  of  experi¬ 
mental  work  (-*-'2, 3, 4, 5)  ^as  £een  perforrae(a  wj.th  metal  foams  or 
loosely  packed  powders  at  high  pressures. 

Theoretical  descriptions  of  the  shock  propagation  in  porous 

solids  are  represented  by  the  work  of  Thouvenin, ^  Butcher ^ 

(1  5) 

and  Kormer,  et  al.  ’  Thouvenin  represents  the  porous 
solid  in  a  one  dimensional  sense  as  a  series  of  plates  and 
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voids  (plate-gap  model) .  The  disturbance  or  shock  is  trans¬ 
mitted  through  the  material  by  shocks  through  the  solid  plates 
and  by  the  material  crossing  the  gaps  at  the  free  surface 

velocity  corresponding  to  the  pressure  state  in  the  solid 
(4) 

plate.  Butcher  has  described  the  porous  solid  m  the  same 
fashion  but  has  restricted  his  application  to  low  stress  regimes. 

(1  j) 

The  Russian  workers,  Kormer,  et  al.,  '  and  Krupnikov, 

(3) 

et  al.,  have  employed  their  very  high  dynamic  pressure 
technology  to  porous  metals  and  have  therewith  constructed  equa¬ 
tions  of  state  for  there  metals  that  describes  off-Hugoniot  bi- 
havior. 

The  basic  approach  taken  by  the  Russian  workers  has  been  ;.o 
assume  an  equation  of  state  of  the  general  form: 

P  a  P  +  P.  +  P  and  (1) 

c  Z  e 

E  =  Ec  +  E£  +  Ee  (2) 

Where  the  subscript  c  denotes  the  0°K  isotherm  and  £  and  e 
are  thermal  components  arising  from  the  lattice  and  electron 
contributions  respectively.  Determination  of  the  0°K  isotherm 
is  described  by  Kormer,  et  al. This  isotherm  Pc  (o)  at  0°K 
is  believed  to  be  valid  over  the  entire  volume  range  from  a 
highly  compressed  solid  in  the  Thomas- v’ermi-Dirac  regime  to 
"cold  ideal  gas"  states.  The  equation 

'.w=IVju1  <3> 

i=l 

is  used  by  references  1,  2,  and  3,  where  the  coefficients,  a^ 
are  determined  by  such  constraints  as  the  heat  of  vaporization. 
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bulk  modulus,  second  and  third  order  elastic  constants,  the 
shock-determined  Gruneisen  ratio,  and  the  Thomas-Fermi-Dirac 
theory.  Sigma,  o,  is  the  state  of  compression  or  p/pQ.  The 
cold  compression  curve  Pc  (a)  together  with  the  lattice 

(a,T)  and  the  electron  P&  (cr,T)  temperature  contributions 
is  used  to  map  states  for  copper  in  P,  V,  E  or  T  thermodynamic 
planes . 

To  provide  the  equation  of  state  with  a  large  enough  range  an 
interpolation  is  required  for  the  lattice  terms.  The  equations: 


\  -  »E  «-*> 


et  =  rn  §■  R  (T“T* 
£> 


cv  -  I  R  [l  +  (l  +  Z)' 


1  -  1  (3y  +  z ) 
4  3  {2  +  Z) 


are  chosen  where 


Z  =  LRT/C' 


C  =  — £ 
c  dp 


is  the  sound  speed  along  the  0°K  isotherm.  L  is  an  experi¬ 
mentally  determined  constant,  R  is  the  universal  gas  constant. 


r~ 


MANUFACTURING  DEVELOPMENT  .  G  £  N  £  R  A  l  MOTOR 


s  CORPORATION 


MSL-7I-30 


where 


CvdT 


and 


(11) 


*o  “  300 °K, 

As  the  parameter  2  -*•  o  ecnaF-s 

Mio  r  ..  .  u  Rations  (4)  and  (5)  reduce 

Mi e-Gruneisen  equation: 


to  the 


P  -  P  e 
Similarly 


c  *  YP  (£“Ec) 


(12) 


°v  3R  a”d  *i  *  V  <*«  Griineisen  ratio. 
The  functional  for*  of  y  used  throughout  is: 

v  .  i  r  »  d2Vd«2 

3  t  ar7d5— 


(13) 


As  Z  oa  * 


p '  »RT- c-  -h  l  i. 


v  2  '  ^  =  3 >  i.e.  the  ideal 


gas  values. 


PTe  ~  9P  ^  In  cosh  (^£) 
2 


r3T, 


and 


=  b 


T  “  ln  c°sh  (|S) 


(14) 


(15) 
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where 


g  = 


din  6 
din  p 


or 


(16) 


6  =  60/o9 


(17) 


The  experimental  results  from  the  porous  work  are  used  to 

establish  values  for  g  and  L  (which  are  assumed  constants) - 

Values  for  6  and  p  are  found  from  measurements  near  absolute 
o  c 

zero  temperatures.  Putting  these  terms  together  then  we  have: 


P 


(p)  + 


3y(p)  +  Z  (P/T) 
1+  Z  (p,T; 


A 

pR (T-T) 


+ 


gp 


b 

MpT 


In  cosh 


#6 (p)  T  % 

<— B - > 


(18) 


E  =  (p) 


,  2  +  Z(P>T)  3 
1  +  Z  (P/T)  2 


R  (T-T)  + 


b‘ 

STpT 


In  cosh  (•— ) 


(19) 


The  procedure  we  used  to  determine  Hugoniots  is  as  follows. 
Equations.  18  and  19  are  not  solvable  for  '£,  so  by  selecting 
a  series  of  temperatures  a  family  of  isotherms  was  calculated 
over  a  range  of  values  from  V  =  0.05  to  0.3  cm3/gm.  The 
pressure  and  specific  energies  for  each  volume  state  along 
each  isotherm  were  plotted.  Next  a  family  of  straight  lines 
representing  E/P  ratios  allowed  by  the  Hugoniot  equation: 


AE  =  E 


Eo  *  1  P  <mVo  ~  V) 


(20) 
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where  drawn  in  the  E  -  P  plane.  Here  m=p/p  ,p  =  porous 
density,  VQ  =  p—  and  V  is  the  shock  state  on  the  Hugoniot. 

In  determing  E/P  only  V  is  allowed  to  vary  ever  the  range  of 
interest  for  each  value  of  m.  Thus  the  intersection  of  E/P 
(at  some  V^)  with  the  line  along  that  same  through  the 
various  isotherms  uniquely  establishes  the  pressure,  energy 
and  temperature  conditions  for  the  Hugoniot  state.  This 
method,  though  graphical,  gives  quite  good  agreement  with  the 
data  reported  by  Kormer  er.  al.  ^ 

A  listing  of  the  computer  routine  used  to  generate  the  required 
isotherms  is  given  in  the  Appendix. 

As  mentioned  previously,  experiments  have  been  performed  at 
high  pressures  with  metal  foams  snd  loosely  packed  powders. 

The  data  from  these  experiments  are  limited  to  small  disten¬ 
tions  and  can  be  described  by  either  the  plate-gap  model  or 
the  Russian  equation  of  state.  The  program  described  here 
was  intended  to  provide  experimental  shock  compression  data 
for  porous  copper  at  high  distentions  up  to  about  0.9  Mb. 

The  results  of  these  experiments  would  then  be  compared  to 
the  predictions  of  the  available  theories. 
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SECTION  II 

DESCRIPTION  OF  EXPERIMENTS 


EXPERIMENTAL  CONSIDERATIONS 

For  this  study  an  "accelerated  reservoir"  light-gas  gun  was 
used  to  launch  flyer  plates  at  velocities  nearing  8  kilo¬ 
meters  per  second.  The  light-gas  gun  range  and  basic  in¬ 
strumentation  have  been  described  in  several  other 
papers  .  This  method  of  experimentation  offers 

significant  advantages  over  the  explosive  techniques  pre¬ 
viously  used.  Unshocked,  stress-free  impactors  of  dis¬ 
similar  material  to  the  specimen  can  be  impacted  over  a 
wide  and  continuous  velocity  range. 

By  varying  the  impact  velocity  and  the  impactor  impedance 
and  measuring  the  shock  propagation  parameters  of  the  target 
a  shock  adiabat  may  be  determined.  This  single  approach  has 
been  employed  to  establish  the  high  pressure  Hugonicts  for 
a  number  of  materials  as  in  Reference  8.  These  shock  wave 
parameters  are  the  snock  velocity  and  the  mass  velocity  or 
particle  velocity  behind  the  shock  front.  In  addition  to 
ewe  va'iiies ,  several"  assumptions  are  implied  in  the 
Rankine-Hugoniot  equations  used  to  describe  the  effect  on 
the  material  of  the  shock  wave.  The  equations  are: 

PH  -  po  -  po  us  “p  (21) 

jr<0  -  1  -  Up/Us  (22) 

*B  -  Ec  m  1/2  <PH  -  Po>  "o  -  V)  !23) 
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where  the  subscripts  o  and  H  refer  to  the  conditions  in  front 
of  and  behind  the  shock  front  respectively  and 

P  is  pressure 

p  =  1/V  or  density 

V  is  specific  volume 

Us  is  shock  velocity 

^  is  particle  velocity 

E  is  specific  internal  energy 

These  equations  assume  a  discontinuity  in  P,  p  and  E  across 
the  shock  front  and  that  PH,  pR  and  Efi  represent  thermodynamic 
equilibrium  conditions. 

The  two  standard  techniques  employed  by  our  laboratory  both 
involve  accurate  and  precise  measurements  of  the  impactor 
velocity.  '  Technique  I  is  used  when  the  target  and  im¬ 
pactor  material  Eire  identical.  In  this  case  u  =  1/2  v.  ,  v.  = 

P  11 

impact  velocity.  We  have  previously  established  that  the 

tially  unshocked  and  thermodynamically  at  the  initial  state. 
Technique  II  is  used  for  dissimilar  material  impacts  and 
depends  upon  detailed  knowledge  of  the  Hugoniot  equation  of 
state  of  the  impactor  material.  The  process  of  establishing 
the  Hugoniot  standards  for  these  experiments  is  Technique  I. 
Even  though  any  given  Hugoniot  state  of  a  standard  is  precise 
to  ^  0.5%,  because  of  the  statistical  nature  of  the  descriptive 
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equation  {usually  a  linear  least  squares  data  fit)  the  re¬ 
sulting  U  (u  )  equation  is  only  accurate  to  about  2%.  The 
®  P 

use,  then,  of  Technique  II  usually  results  xn  Hugoniot  data 
ill  the  unknown  material  of  an  accuracy  of  ^  2-4%.  The  as¬ 
sumption  of  instantaneous  state  discontinuity  (ignoring  low 
pressure  elastic  waves)  ensures  that  the  high  pressure  shock 
process  is  truly  adiabatic,  i.e.,  all  the  change  in  internal 
energy  occurs  from  compression  alone,  and,  without  extremely 
difficult  stepwise  analysis,  is  the  only  way  to  describe  the 
shock  on  a  macroscopic  scale.  In  addition,  the  assumption 
of  thermodynamic  equilibrium  establishes  the  only  basis  at 
present  for  analytical  treatment  of  the  shock  data  from  end- 
pojlnt  to  endpoint.  If  either  of  these  assumptions  fail  at 
high  pressures,  the  Hugoniot  interpretation  of  the  measurements 
also  fails.  There  are  several  possible  ways  this  might  be  ob¬ 
served. 

a.  Konsteady  shock  prcpagation.  Since  it  is  obvious  that 
shock  waver  do  attenuate  the  criterion  must  be  established 
at  some  limit.  I rimarixy  very  high  pressure  shocks  atten¬ 
uate  geometrically,  but  nany  materials  evidence  irreversible 
energy  losses  through  pnase  transitions  both  solid/solid 
type  and  solid/liquid  or  solid/liquid/gas.  This  is  ob¬ 
servable  by  non-constant  shock  propagation  velocities  with 
increasing  specimen  thickness. 

b.  Cross  correlation  of  state  measurements.  Shocks  might 
exhibit  non-equilibrium  conditions  if  the  experimental 
redundancy  of  recording  systems  are  accurate  enough.  An 
example  right  be  the  simultaneous  measurement  of  the  shock 
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pressure-time  history  and  shock  velocity,  or  some  other 
shock  parameter  which  would  allow  redundant  calculation 
of  the  shock  state. 

c.  The  appearance  of  structured  or  multiple  waves .  Multiple 
waves  at  very  high  pressures  may  indicate  the  existence 
of  non-equilibrium  conditions  behind  the  shock,  assuming 
that  phase  transitions  are  accounted  for. 


MATERIAL  DESCRIPTION 

The  distended  copper  employed  as  test  material  was  produced 
by  a  process  developed  at  GM  Research,  General  Motors  Technical 
Center,  Warren,  Michigan.  Known  as  MetNet,  the  specimens  bulk 
properties  are  shown  in  Table  I.  Microscopic  examination 
reveals  chat  the  copper  MetNet  structure  is  web like.  The 
webbing  is  the  interstitial  volume  between  tightly  packed 
spheres  of  nearly  identical  diameter  as  shown  in  Figure  1. 

Each  cell  interconnects  with  its  neighbors  through  nearly  cir¬ 
cular  openings  -  the  number  of  openings  depending  on  the  dia¬ 
meters  of  the  cell  and  its  neighbors.  Openings  per  cell  range 
from  10  to  15.  The  cross  section  of  the  intracellular  web  is 
triangular  with  the  three  side*,  concave  outward,  and  a  similar 
triangular  hollow  center.  The  geometry  remains  identical  with 
increasing  density  only  the  cross  sectional  thickness  of  web 
increases . 

The  target  densities  of  approximately  33%,  22%  and  11%  of  solid 
copper  were  selected  as  compromises  between  the  constraints  of 
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TABLE  I 

PHYSICAL  PROPERTIES  OF  COPPER  METNET 


Sanqole  Average  Density 


3 

(gm/cm  ) : 

Distention: 

2.94 

3.08 
rms  0.03 

1.99 

4.46 
rms  0.03 

1.001 

8.74 
rms  0.08 

* 

Spectro/Chemical  Analysis 

Element  %_Comp. 

Copper 

95.02 

98.30 

96.99 

Nickel 

0.92 

1.49 

2.38 

Zinc 

Trace 

Trace 

0.43 

Iron 

Trace 

0.10 

Trace 

Pore  Size 

Range  (in. ) : 

0.008"— 

0.035" 

0.010"- 

0.040" 

0.010 

0.040 

Approx.  Ave.  (in.) 

0.012" 

0.015" 

0.015 

Sound  Speed 

Approx. 

Longitudinal  Vel. 

1.92 

1.55 

0.98 

(mm/vtsec) 


Tests  perforated  by  Met-Chem  Testing  Laboratories,  Inc. 
20062  John  R,  Detroit,  Michigan. 
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the  manufacturing  process  —  i.e.  limited  maximum  densities  of 
•u  50%  solid,  and  mechanical  requirements  for  handling.  MetNet 
can  be  produced  at  distentions  of  100  or  more. 

The  specimen  densities  and  porosities  were  determined  by  two 
methods:  a)  water  saturated,  dry  and  immersed  weights  and 
b)  bulk  dimensions  with  dry  weights. 

All  target  and  test  specimens  were  taken  from  single  10"  x  12" 
sheets  1/4"  thick  for  each  distention.  Density  measurements 
were  made  on  each  target  specimen  as  well  as  on  several 
selected  samples  from  locations  throughout  the  respective 
sheets.  The  water  immersion  method  was  employed  and  by  soaking 
the  specimens  overnight  in  a  pressurized  water  bath,  the  bulk 
density  and  porosity  measurements  were  within  0.5%  of  the 
values  determined  geometrically.  Individual  bulk  density 
measurements  may  tend  to  be  in  error  due  to  the  small  sizes 
for  some  of  the  target  specimens,  nevertheless,  the  measured 
values  shown  in  Table  II,  Til,  and  IV  are  employed  in  the 
Hugoniot  calculations  - 

The  nickel  contaminant  shown  in  Table  I  is  a  result  of  the 
manufacturing  process  for  MetNet  but  is  sufficiently  low  in 
view  of  the  great  distentions  to  not  be  a  significant  factor 
in  the  high  pressure  behavior. 


SPECIMEN  PREPARATION 

Target  specimens  were  core  drilled  from  the  stock  sheets 
slightly  oversize.  Turning  to  final  test  diameter  was  per¬ 
formed  after  one  face  of  the  target  was  machined  flat.  It 
was  learned  that  although  the  machining  operations  tended  to 
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smear  and  distort  the  web  network,  the  deformed  zones  could  be 
removed  by  lapping  since  they  were  only  about  three  pore  dia¬ 
meters  thick.  The  lapped  surfaces  appeared  microscopically 
to  be  nearly  ideal  cross  sections  of  the  foam.  For  the  most 
part,  the  target  edges  were  left  in  the  smeared  condition 
since  this  enhanced  target  mounting  by  epo::v  casting.  Care 
was  taken  to  seal  the  specimen  edges  so  that  the  mounting  epoxy 
did  not  seep  into  the  target  thus  changing  its  apparent  bulk 
density. 

The  planarity  of  the  lapped  target  faces  was  checked  on  a 
Zeiss  light  section  microscope.  While  the  surface  roughness 
was  of  the  order  of  the  cellular  diameters,  the  sectioned 
webbing  was  uniform  and  the  impact  and  rear  surfaces  coplanar 
to  within  ^  0.0001  inch. 

Oxidation  and  corrosion  of  the  copper  was  kept  to  a  minimum 
by  storing  the  specimens  in  a  vacuum  desiccator  after  they 
were  annealed  in  a  reducing  atmosphere. 

Target  handling  was  a  problem  only  for  the  highest  distention 
MetNet  which  could  be  deformed  by  finger  tip  pressure.  Jlount- 
ing  and  machining  fixtures  and  rareful  handling  prevented  in- 
adver'iant  marring  collapsing  of  the  test  specimens. 


TARGET  DESIGN 


Pin  Targets 


The  majority  of  the  targets  used  in  these  experiments  were 
of  the  four  pin  design  described  in  earlier  papers .  The 
instrumentation  of  these  targets  permitted  the  measurement  of 
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two  shock  transit  time  intervals  during  each  shot,  as  well 
as  a  recording  of  the  impactor  tilt  relative  to  the  pin  axis. 

Several  5  cm  diameter  targets  were  instrumented  wich  mica 
sheathed  ipanganir;  wire  piezoresi stive  transducers^^.  Front  and 
rear  surface  pins  were  also  used  to  measure  the  shock  transit 
time.  These  targets  consisted  of  two  disks  t,  4.5  m  thick 
with  the  sheathed  gage  sandwiched  in  between.  Although  several 
versions  of  sheathing  and  assembly  were  constructed,  none  of 
them  successfully  recorded  a  pressure  profile. 


Optical  Targets 

Several  targets  were  designed  to  be  optically  instrumented  by 
a  Beckman  and  Whitely  Model  339B  high  speed  continuous-writing 
streak  camera.  Three  wedge  shaped  targets  were  tested  at  the 
three  distentions  for  nonuniform  shock  propagation  with  thick¬ 
ness  at  relatively  low  pressure.  Target  design  is  shown  in 
Figure  2.  Light  reflecting  from  the  front  surface  mirror/ 
target  interface  is  extinguished  upon  the  arrival  of  a  shock, 
thus  photographically  permitting  continuous  shoe,  transit 
time  interval  determination.  A  record  obtained  from  these  shots  is  shown 
in  Figure  3  graphical  reductions  of  Shots  401,  407  and  408  are  shown  in 
Figures  4,  5,  6  and  7. 

Single  and  double  step  targets  were  constructed  for  use  at 
higher  impact  velocities  than  the  three  shots  above.  The 
design  for  these  targets  is  shown  in  Figure  8.  Over  the  range 
of  impact  velocities  achievable  by  the  29  mm  launch  tube  all 
the  inpact  events  are  self-luminous,  i.e.,  the  impactor 
strike  at  the  target  front  surface,  and  the  arrival  of  the 
shock  at  the  rear  surface  glass  witness  plate  are  light  emit¬ 
ting  events.  Figure  9  displays  one  such  experiment  record. 


IS 


WEDGE  SHOT  NO-  401 


SHOCK  VELOCITY  MM/USEC 
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Figure  6  Optical  Wedge  Target,  Shock  Velocity 
vs.  Target  Thickness,  Shot  No.  407 


SHOCK  VELOCITY  MM/USEC 


TARGET  THICKh£5S 


Figure  7  Optical  Wedge  Target,  Shock  Velocity 
vs.  Target  Thickness,  Shot  No.  40S 
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GLASS  WITNESS 
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EPOXY  SEALANT 


PLEXIGLAS  RING 
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Figure  8  Design  of  Optical  Double  Stepped 

Target  for  High  Pressure  Experiments 
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A  single  optical  test  was  performed  which  was  also  instrumented 
by  pins.  This  shot  was  designed  to  cross  check  the  correlation 
with  the  pin-only  target  data  to  the  optics-only  target  results. 


DATA  ANALYSIS 

The  basic  parameters  measured  in  these  experiments  were  the 
impact  velocity  of  the  impactor  and  shock  velocity. 

The  impactor  velocity  measuring  system  consists  of  a  laser 

(8) 

triggering  system  and  two  short  duration  flash  X-rays. 

With  this  system,  impactor  velocities  are  measured  accurately 
to  0.05%.  The  triggering  system  consists  of  a ■ helium-neon 
gas  laser  aimed  at  a  photo-detector  across  the  impact  chamber, 
orthogonal  to  and  intersecting  the  line  of  flight  of  the  pro¬ 
jectile.  A  photomultiplier  monitors  the  laser  light  output 
through  a  set  of  masks  and  a  narrow  bandpass  optical  filter. 
When  light  interruption  occurs  due  to  projectile  passage,  a 
sharp  change  of  voltage  level  is  converted  into  a  signal  of 
sufficient  amplitude  to  trigger  a  Field  Emission  Corporation 
30  nanosecond  dual  flash  X-ray  unit.  The  X-ray  flash  exposes 
a  Polaroid  film  plate  on  the  opposite  side  of  the  chamber  by 
means  of  a  fluorescing  intensifier  screen.  The  trigger  and 
X-ray  flash  system  is  then  duplicated. to  record  the  passage 
of  the  projectile  in  the  second  field  of  view  30.5  cm  further 
down  range. 

The  spacing  between  the  two  X-ray  field  centerlines  is  indica¬ 
ted  by  fiducial  wires  which  are  measured  by  an  optical  com¬ 
parator  to  within  0.2  mm.  Measurements  of  the  impactor  face 
position  relative  to  the  window  fiducials  allow  calculation 
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of  actual  projectile  position  and  travel  over  the  time  interval 
measured  between  flash  exposures. 

Measurement  of  shock  velocities  was  accomplished  through 
recording  of  the  time  interval  between  impact  and  shock  arrival. 
The  closure  of  coaxial  self-shorting  pins  provided  signals  to 
start  and  stop  Eldorado  Model  793  time  interval  meters.  These 
counters  have  a  specified  time  resolution  of  +  1/2  nsec  and 
are  read  digitally  to  the  nearest  nsec. 

The  optical  records  were  converted  to  time  intervals  through 
reduction  from  analog  to  digital  form  using  an  optical  com- 
paritor.  The  digital  records  are  then  computer  analyzed  into 
time  intervals  and  shock  wave  velocities.  Corrections  for 
impactor  tilt,  target  surface  non-parallelism,  and  statistical 
averaging  are  employed  in  the  coiqputer  analysis.  The  actual 
measurements  are  relative  time  positions  on  stationary  film 
in  the  continuous  writing  streak  camera.  These  relative  film 
event  positions  are  converted  to  time  intervals  in  terms  of 
the  camera  writing  rate  through  knowledge  of  the  turbine  an¬ 
gular  velocity  and  corrections  are  applied  for  the  optical 
and  mechanical  nonlinearities. 

Although  all  the  data  here  are  presented  as  shock  velocities, 
it  should  be  kept  in  mind  that  the  actual  recorded  data  is  in 
the  form  of  relative  event  times  and  velocities  are  inferred 
on  the  assumption  of  constant  velocities. 
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SECTION  III 

EXPERIMENTAL  RESULTS 


This  section  presents  the  results  of  the  data  analysis  pre¬ 
viously  described.  A  few  prelirainary  remarks  are  necessary. 
The  attests  to  record  the  press  are- time  profiles  of  the 
three  distentions  of  MetNet  were  unsuccessful.  All  of  the 
sheathed  manganin  r»ire  transducers  imbedded  in  the  targets 
shorted  immediately  upon  the  arrival  of  the  transmitted  wave. 
We  have  no  experimental  evidence  en:bling  us  to  ascertain 
whether  the  shorting  occurred  along  the  gage  leads  or  in  the 
active  area.  Even  additional  buffering  by  increasing  the 
sheath  thickness  had  no  apparent  effect  on  gage  lifetime. 
Interpretation  of  the  three  wedge  -targets  at  the  relatively 
low  impact  velocity  of  ^  2.85  loVsec  was  that  no  appreciable 
wave  propagation  velocity  change  occurs  in  the  thickness 
range  fror.  2  to  10  mm.  Thus,  a  single  thickness  of  i  5  nan 
was  employed  for  most  of  the  pin  experiments. 


PIN  DATA 

Compilations  of  the  pin  experimental  results  are  presented 
in  Tables  II,  III  and  TV  for  copper  MetNet  at  the  average 
distentions  of  mave  =  3.08,  4.46  and  8.74  respectively.  These 
data  are  displayed  simultaneously  in  the  shock  velocity-particle 
velocity  plot.  Figure  10.  Simple  curves  drawn  through  the 
data  suggests  that  experimental  scatter  is  quite  significant. 
Interpreting  the  scatter  as  curvative  is  rather  dubious  since 
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conversion  to  the  pressure-volume  plane  results  in  double  and 
possibly  triple  valued  states  at  some  values  of  V.  This  may 
be  seen  in  Figure  15  by  simply  connecting  the  dots. 


OPTICAL  DATA 

Tables  V,  VI  and  VII  present  the  numerical  results  from  the 
optical  records  as  well  as  all  the  pin  data.  A  different 
'nterpretation  had  to  be  made  of  these  records  than  for  the 
pin  records.  The  recording  time  of  these  optical  experiments 
was  much  longer  than  with  the  pin  technique.  That  is,  optical 
recording  of  events  is  possible  as  long  as  the  glass  witness 
plate  remains  undestroyed.  These  glass  witness  plates  were 
not  silvered  so  that  events  prior  to  the  actual  rear  surface 
arrival  were  observable.  Figure  11  is  a  record  of  Shot  521 
into  a  target  with  no  step  and  *  9  mm  thick,  showing,  1)  the 
arrival  at  the  target  front  face,  2)  luminousity  from  inside 
the  open  celled  structure,  3)  arrival  at  the  witness  plate  of 
the  luminous  streams,  4)  arrival  at  the  interface  of  a  des¬ 
tructive  wave,  and  5)  arrival  of  that  destructive  wave  at  the 
rear  surface  of  the  witness  plete. 

Events  3  and  4  are  characterized  more  as  zones  of  changing 
luminousity.  An  indication  of  the  substantial  energy  represented 
by  the  discrete  luminous  streams  is  the  shock  diamond  pattern 
in  the  glass  witness  plate  just  after  Event  4.  The  time 
interval  of  Event  4  is  identifiable  in  all  the  optical  records 
and  is  significantly  larger  than  that  recorded  by  the  self¬ 
shorting  pins.  In  fact,  the  pin  recorded  time  intervals  for 
nearly  identical  inpact  conditions  (and  for  Shot  527  identical 
inpact  conditions)  corresponds  to  the  optical  Event  3.  In 
addition,  the  scatter  ia  the  pin  records  correlates  well  with 
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the  optically  observed  irregularity  of  stream  arrivals. 

Tables  V,  VI,  and  VII  are  also  conpilations  of  all  the  experi¬ 
mental  work  at  the  three  MetNet  distentions.  The  shock  wave 
velocity  data  are  correlated  to  copper  impactor  velocity  and 
segregated  by  recording  technique  and  recording  station  posi¬ 
tion.  The  impact  velocities  contained  in  parentheses  are 
confuted  values,  i.e.  the  velocity  a  copper  impactor  would 
have  to  induce  the  same  state  in  the  target  that  was  achieved 
by  the  actual  Fansteel-77  impact.  The  optical  data  columns 
contain  the  recorded  velocities  of  both  the  first  and  the  final 
"waves",  and  in  the  cases  of  the  stepped  targets  both  thick¬ 
ness  station  values  are  listed  separately  but  under  the  same 
shot  number.  The  approximate  time  interval  between  the  two 
waves  is  listed  in  column  "Average  At."  The  pin  data  are 
tabulated  under  either  station  A  or  B  and  the  station  thick¬ 
nesses  are  listed. 

Figures  12,  13  and  14  are  shock  velocity- impact  velocity  plots 
of  the  data  contained  in  Tables  V,  VI,  and  VII  respectively. 

Shot  numbers  are  included  to  identify  the  impact  velocities 
for  the  experiments.  The  endpoints  of  the  error  bars  show  the 
first  and  last  wave  velocities  for  the  optical  shots  and  the 
data  scatter  for  the  pin -shots. 
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SECTION  IV 

DISCUSSION  OF  RESULTS 


COMPARISON  OF  PIN  DATA  TO  RUSSIAN  WORK 

A  primary  objective  of  the  program  described  here  was  to 
cciqpare  experimentally  similar  materials  testing  between  our 
efforts  at  large  distentions  and  the  Russian  Work^'"^  with 
distended  copper  at  m  =  1.57/  2.0,  3.01  and  4.0.  Our  disten¬ 
tions  were  m=  3.08,  4.46  and  8.74.  The  experiments  at 
m  =  3.01  and  3.08  were  expected  to  overlap.  Since  essentially 
identical  experimental  techniques  were  employed  (except  for 
irapactor  velocities),  i.e.  electronic  pin  switches,  similar 
target  size,  time  interval  measuring,  any  significant  dif¬ 
ferences  in  results  can  only  be  attributable  to  the  specimen 
structure  itself.  As  can  be  seen  in  Figures  15  and  16,  there 
is  considerable  difference  between  the  pin  recorded  states  of 
this  and  the  Russian  work.  Dashed  lines  are  used  to  indicate 
the  predicted  Hugoniots  for  the  values  of  m  shown.  The  pin 
data  plotted  for  these  distent,  .ons  also  deviate  significantly 
from  the  predicted  Hugoniots.  The  Russian  experimental  data 
are  plotted  with  open  symbols  as  shown  in  the  legend. 


QUALITATIVE  DISCUSSION  OF  OPTICAL  RECORDS 

The  additional  observations  performed  by  us,  i.e.  the  optical 
tests,  suggest  a  mechanism  for  the  above  discrepancies.  From 
the  fact  that  our  records  show  a  displacement  in  the  P-V  plane 
to  the  right  of  the  predicted  Hugoniots  it  can  be  inferred  that 
our  recorded  shock  velocities  are  too  high.  Indeed,  the  optical 
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records  suggest  this  also.  Event  #3  in  Figure  11  has  been 
interpreted  as  the  first  arrival  of  a  pressure  pulse  both 
because  it  represents  a  zone  of  extinguished  light  and  because 
its  time  interval  corresponds  closely  to  that  recorded  by 
self-shorting  pins.  An  additional  implication  suggested  by 
the  emitted  light  from  Event  #2  is  that  Event  #3  is  the  arrival 
of  a  hot  high  velocity  vapor  or  molten  stream.  The  open  cell 
structure  of  MetNet  would  probably  allow  fine  streams  to 
propagate  several  cell  diameters  before  encountering  further 
stationary  solid  material.  The  effect/  however/  of  using  the 
final  wave  or  Event  #4  as  the  stable  equilibrated  shock  is 
to  shift  the  P-V  states  so  fair  that  some  states  lie  to  the  left 
of  solid  Hugoniot.  Some  precompression,  then,  must  be  occurr¬ 
ing  in  the  first  wave,  at  least  to  the  extent  of  modifying  the 
initial  specific  volume  an  unknown  and  indeterminable  amount. 

The  optical  records  then  invalidate  the  Hugoniot  states  deter¬ 
mined  by  the  simplistic  pin  experiments  by  destroying  the 
assumption  of  shock  state  equilibrium  and  through  lack  of  know¬ 
ledge  of  the  intermediate  and  final  volume  states. 

In  Figure  17  the  time  vs  distance  data  from  the  optical  records 
of  step  target  inpacts  are  presented.  For  all  three  distentions, 
the  first  wave  is  seen  to  parallel  the  second  wave  between  the 
steps.  Also,  the  waves  appear  to  be  decreasing  in  velocity 
with  distance  traveled.  It  is  not  possible  to  determine  if 
this  is  due  to  an  initial  instability  or  if  the  waves  are 
attenuating  with  distance.  Finally,  for  the  three  distentions 
shown,  the  wave  trajectories  are  similar  for  similar  impact 
velocities. 

It  appears  that  the  first  wave  and  second  wave,  after  estab¬ 
lishing  some  separation,  propagate  at  the  same  velocity.  The 
separation  of  the  two  waves  we  interpret  to  be  associated  pri¬ 
marily  with  the  cell  size  and,  to  some  minor  extent,  the  dis¬ 
tentions  . 
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If  the  waves  are  decreasing  in  velocity  with  distance, 
calculation  of  states  at  a  point  cannot  be  made  solely  on  the 
knowledge  of  transit  times. 

T.ie  wave  trajectories  for  the  three  distentions  are  similar 
for  similar  impact  velocities  as  seen  here  and  also  as  inferred 
from  the  wave  velocities  vs  inpact  velocity  as  displayed  in 
Figures  12,  13  and  14. 

No  other  experimental  techniques  at  these  high  pressure  con¬ 
ditions  have  been  reported  which  could  have  detected  the 
existence  of  a  double  wave.  This  is  particularly  true  if  the 
leading  wave  carries  sufficient  energy  as  suggested  in  these 
experiments  to  significantly  alter  the  stationary  materials 
state.  That  this  leading  wave  phenomenon  is  peculiar  to 
materials  such  as  MetNet  is  a  moot  point.  It  seems  difficult 
to  ignore  the  possibility  of  a  ramped  non-equilibrium  wave 
even  in  finely  powdered  porous  specimens  with  any  significant 
distention. 


44 


ANUrACTURING  DEVELOPMENT  •  GENERAL  MOTORS  CORPORATION 


MSL-71-30 


SECTION  V 

CONCLUSIONS 


On  the  basis  of  the  experimental  results  described  above,  the 
following  conclusions  have  been  reached: 


1.  Optical  observation  of  large  distention,  open  celled  cop¬ 
per,  showed  a  multi-wave  complex  loading  structure  at  high 
pressures . 


2. 


3. 


There  is  some  evidence  that  the  waves  in  the  distended 
copper  are  slowing  down  with  distance  of  propagation. 
This  is  based  on  the  data  presented  in  Figures  6,  7 
and  17.  The  wedge  data  showed  variations  in  velocity 
of  1  to  3%  at  low  pressures.  The  optical  step  data 
indicated  a  change  in  velocity  between  steps.  This 
may  be  an  initial  instability  or  a  continuous  reduction 
of  velocity  with  distance  of  propagation. 
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To  a  first  approximation,  the  same  wave  speeds  were  ob¬ 
served  in  samples  of  different  distentions  for  the  same 
impact  velocities. 


4. 


Comparison  of  the  results  of  the  experiments  with  the  pre¬ 
dictions  of  the  Russian  equation  of  state  is  poor  due  to 
the  multi-wave  structure  we  observed. 
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In  summary,  the  analysis  of  the  experiments  with  porous  copper 
show  that  the  loading  process  cannot  be  represented  as  a  sinrole 
shock  process  even  at  the  highest  pressures.  A  description 
of  the  process  would  have  to  include  the  cell  size  in  addition 
to  the  porosity. 
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**  POROUS  MATERIALS  EQUATION  OF  STATE  PROGRAM 

'•  DIMENSION  V(50) ,S( 50) ,PCOLD< 50) .ECOLDt 50) ,PTL( 50)  ♦ETLISOl * 
1P1  E(  50)  *  ETF  (.50  )  .RETA(5C)»Z(50)»CCK(50)  »GAM(50)  ,  P(50)  * 

2ECAL  (  50  )  »TMATLf  5)  .HC0S(50)  ,T  (50)  .EnUGt  50)  .CM (50)  »DPDV(5P‘ 
COMMON  A1,A2»A3*A4»A5»A6»A7 


INPUT  CONSTANTS 

0  DEG.  <  SOLID  DENSITY  - 
DISTENTION  - 
SPECIFIC  HEAT  (V)  - 
SPECIFIC  HEAT  <P)  - 
INITIAL  TEMPERATURE  - 
ELECTRON  SPECIFIC  HEAT  - 
ELECTRON  TER"  - 
7  MATERIAL  CONSTANTS  - 
START  VOLUME  - 
FINAL  VOLUME  - 


RHOSO 

DM 

CV 

CP 

TO 

BETAO 

G 

A1  THRU  A 7 

VST 

VFI  ii 


RFAD(  5»10)  (TMATLdJ  .1  =  1,5) 

READ (5,11) RHOSO. DM, CV.CP.Tll) .BETO.G 
READ(5,11)A1,A?,A3.A4.A5.A6.A7 
RE  AD ( 5 , 1 2 ) VST , VF I N 

10  FORMAT (5A4) 

11  FORMAT (7F1 0.5) 

12  FORMAT  (  2F10.5  )  _ _ _ _ 

tTSSSSTm*®’ 

DO  40  K  =  1  ♦  5  0  ava»)ab)e_c_H7_^^ 

R=8.3l432E+07/63.54 
D=12.CE+06 
TTK=10C. - 
VSO=l. /RHOSO 
DELV= (VFIN-V5T) /50. 

V .  1 ) =VST 

Si  1  )=VS0/V( 1 ) 

DO  100  1=2.50 
V  (  I ) =V ( 1-1 ) +PELV 
S<  I  )=VS0/V( I ) 

CALL  ISOTHiVSO.SU  )  ,PCOLD(  I  )  .ECULDU  )  »GAM(  I  )  .DPDVl  I  )  ) 
PCCLD ( I ) =PCCLD( I ) *1 • OE+12 
ECOLD ( I ) -CCOLD1 I>«1.0E^12 
DPDV ( I )=1.0F+i?*DPDV( I ) 

BETA  (  I  )  =RET0  /(Sm**G) 

CCKC I )=  VS0*DPDVI I ) 

I F { V { I )-. 23)20.20. 21 
?0  Z  (  I  )=L*R»T  (.<)  /CC<(  I  ) 

PTL ( I )  =  ( 3.*GAV( 1 )+Z ! I ) )*R*CT (<  J-TTO / 5  ( l.+Zt II)*V(II) 
ETL ( I ) = ( 2.+Z ( I ) )*3.*R*(T(K)~TTK)/( (l.+Z ( I ) )*2.) 

GO  TO  22 

21  DTL(  I  ) =R*{ T  (Kl-TTK.)  /  V(  I  ) 


Preceding  page  blank 
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ETLU  )=R*(T(K)-TTK)*3.0/2. 

22  CONTINUE 

NCOS (I )=COSH(RFTA(I 1 ,T(K) ,B> 

PTE ( I  5  =G*R**2 / ( BETA ( 1 )*V( I  1 )*ALOG(HCOS(  I  1  1 
ETC ( I )  =  ( 8**2/BFTA( I )  I*ALOG(HCOS( I  1  ) 

P ( I ) =PCOLD ( I ) +PTL ( I ) +PTE ( I  ) 

ECAL ( I ) =ECOLD ( I l+ETL ( I )+ET£( I ) 

EHUG( I )=.5*P( I )*(DM*VS0-V( I ) ) 

CV,  ( I  )s(2,*ECAL(I)/P(IJ-VI  I  1 l/VSO 
PI  I )=P(I )*1.0E-06 
°CGLD( I )=PCOLD( I )*1.0E-06 
PTLI I )=PTL( I )*1.0E-06 
pTE  (  I  )=PTE( I )*1.0E-06 
IV'  CONTINUE 

WRITE (6*13) I TVATL ( I  )  .1  =  1,5) 

WR I TE ( 6  *  16 1 RHOSO . OX  * CV .CP . T ( K 1  ,  BETO .  G 
WRITE (6*17) 

00  111  1=1. SO 

111  WRITE  ( 6.14 )  VCD  *SU>  »PU) .ECAL 1 1 ) ,PTl(  I ) »ETL( 1 > »PTE( I ) »ETE  t I  ) 

•c  F0RMATI1H  4F15.6) 

WRITE (6.71 (VII) ,P( I J .£CAL( 1 1 .EHUGI I ) .  Z( 1 1 .1  =  1.501 
7  FORvAT ( lH  .5E15.6) 

T(<+1 )=Tf<)+10000.0 
CALL  DATSw(O.'.SWO) 

IF  ( I SWO-1 ) 40  .200.40 
20 0  UR  I TF (6.9) 

°  F0RVAK1H  »2(« 12HSPEC.  VOLUME  »  5X  *  1 OrtSOU.MD  VEL.  .  IX  .  14HINTERP.  PA.,... 
1 . ,  / / , 5X , PHCM**3/GM . 9X , 6HCM/SEC  » 1 1 X 1H2 . / / 1 

DO  101  1=1.50 
CCK ( I  1 =SCRT ( CCK ( 1 1 1 

.0]  WRITF ( 6 » 19 ) V( 1 1 .CCK ( I  1.2(1) «DPDV( I  1 

WRITE (6.15) (VI  I ) »St I  1 »PCOLD( I  1 .ECOLD( I  1 ,GAM( 11.1  =  1,50) 

13  FORMAT ( 1H1 . 3PHEQUAT ION  OF  STATE  CALCULATIONS  FOR  -  .5A4.//1 

14  FORMAT ( 1H  *  8F1 5  «  6 1 

15  cORvAT { 1H1 .23HISOTHERMAL  STATE  VALUES . //2X , 12HSPEC.  VOLUME. 4X, 

llll-;COvPRESSICN.3X,13H  PRESSURE. 4X ,11H  ENERGY .// ,6X , 

?flHCM**3/GM. 10X .4HV0/V.7X ,8H  R »RS » 8X  .7HERGS/GM » / / ( 5E15. 6 > 1 

16  FORMAT ( 1H0.19H0  DEG.  K  DENSITY  =  »F1C • 5 . 5X « 13HDI STENT  ION  =  »F7.*, 
15X.?0HSP£CIFIC  HEATS  V  = .E15.6 »2X, 3HP  *,E15.6.// 

P  5Sf  15HINITI  ai_  TEv.s.  s.FIO.1  .5X.25HINIT.  ELEC.  SPEC.  HEAT  =  .  E 1 0 . 
3 ' ♦ 5X . 1 7HELEC . RON  COEF.  =  .F10.6.//1 

17  FORMAT (2X» 

1  i 2HSPCC*  vGLU‘-'c  .  4X .  1  IHCG-’.PREJSION.  3X  .  13HTCT  .  PRESSUR,  .-ll«TjT 
2.  ENERGY,  2X«13nLA'.  PRESSURE. 4X.11HLAT-  tNERGY  »  IX  ,  14HELEC  .  .■’■vESSv.v 
3F.3X, 

412FELEC.  FNFRGY  »//6X»fiHCv.**3/GM,13X,4HV0/V»7X,eH  BARS  ,  8X  .  7- 
5S/GM .11X.4HBARS.SX. 7HERGS/GM , 11X .4HBARS , 8X .7HERGS/GM. / / 1 
40  CCNTINUF 
20!  CALL  EXIT 
END 
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